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Abstract 

Motivated by recent anomalous CDF data on Wjj events, we study a possible explanation within 
the framework of the two-Higgs doublet model. We find that a charged Higgs boson of mass ~ 
140 GeV with appropriate couplings can account for the observed excess. In addition, we consider 
the flavor-changing neutral current effects induced at loop level by the charged Higgs boson on 
the B meson system to further constrain the model. Our study shows that the like-sign charge 
asymmetry can be of O(10~^) in this scenario. 
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Recently the CDF Collaboration reported data indicating an excess of Wjj events where 
W decayed leptonically l|. The excess shows up as a broad bump between about 120 and 
160 GeV in the distribution of dijet invariant mass Mjj. This dijet peak can be attributed 
to a resonance of mass in that range, and the estimated production cross section times the 
dijet branching ratio is about 4 pb. However, no statistically significant deviation from the 
standard model (SM) background is found for Zjj events. Events with b-jets in the excess 
region have been checked to be consistent with background. Moreover, the distribution of the 
invariant mass of the iujj system in the Mjj range of 120 to 160 GeV has been examined and 
indicates no evidence of a resonance or quasi-resonant behavior. The D0 Collaboration also 
performed a similar analysis, but found no excess Wjj events ^ . While waiting for further 
confirmation from the Large Hadron Collider at CERN for the result of either experiment, 
it is nevertheless worth pursuing the cause of the anomaly observed by CDF. 

Many papers have discussed different possible explanations for the excess 

IJQ- Most of 

them try to explain the excess by introducing one or more additional new physics particles. 
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. Some others 



3l| , and neutral color-singlet vector boson 
analyze the anomaly considerin g sc alar bosons, such as technipion |^| (including technirho). 



sup er-partners of fermions 
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13[] (fermions in SUSY model are also considered in Refs. [3|, 
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radion |26l | , sc alar 



32| . scalars with flavor symmetry 
doublet with no vacuum expectation value (VEV) 33|, and new Higgs bosons [2 71 . |29| . In 
Ref. 29|, the two-Higgs doublet model (THDM) is discussed with flavor-changing neutral 
current (FCNC) interactions allowed through neutral Higgs boson [H^ and A^) exchanges. 
Their result favors a light charged Higgs boson. However, allowing large Yukawa couplings 
to leptons in their work has the problem that lepton pairs will be copiously produced, which 
is not the case in the CDF data. There are also attempts to explain this puzzle within SM 
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In this letter, we explore another scenario in the THDM as an explanation. The fact that 
the excess dijets are non-b-jets suggests that the new resonance may not couple universally 
to quarks. A scalar particle can accommodate this feature more easily than a gauge particle. 
We show in Fig. [1] two processes in the THDM that can possibly contribute to the excess 
events. The dijets come from the decay of the charged Higgs boson H^. Since the CDF 
Collaboration does not observe any resonance in the invariant mass spectrum of iujj for the 
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excess events, we require that the mass of the pseudoscalar Higgs boson A*^ to be sufficiently 
high. In this case, only Fig. [Tt^a) is dominant, with the mass of the charged Higgs boson 
~ 140 GeV, as suggested by data. Moreover, we assume that the width of is 
sufficiently small in comparison with the jet energy resolution of the experiment. We note 
that this is only one possible scenario in the model. Another scenario is that and are 
interchanged in Fig. [H and so are their masses. We also note in passing that the assumed 
mass of ~ 140 GeV for or Ao is consistent with the lower bounds of 76.6 GeV for 
42 1 and 65 GeV for A^ {41 1 from LEP experiments. 
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FIG. 1: Diagrams contributing to the Wjj events in the two-Higgs doublet model. 



The Yukawa sector of the THDM is given by 



D i 



(1) 



where Hi^2 are two Higgs doublet fields, Hk = i'T2Hl, Ql represents left-handed quark 
doublets, Ur and Dr are respectively right-handed up-type and down-type quarks, and 
Y^2^ are Yukawa couplings. Here we have suppressed the generation indices. The fields 
Hi and H2 can be rotated so that only one of the two Higgs doublets develops a VEV. 
Accordingly, the new doublets are expressed by 



h = sin f3Hi + cos f3H2 



H = cos /3Hi — sin PH2 



G+ 

(i; + /i" + zG0)/V2 



(2) 



where sin/? = vi/v, cos/3 = V2/V, v = ^/vf + fg, (H) = 0, {h) = v/\/2. In our scenario, we 
assume that has mass ~ 140 GeV and is responsible for the excess Wjj events observed 
by CDF. As a result, Eq. ([1]) can be rewritten as 



-C 



with 



Y = QLY^'Unh + QLY''DRh 
+ QlY^'UrH + QlY^DrH 

= sin l3Yf + cos I3Y^^ , 
Y^ = cos l3Yf - sin f3Yl , 



(3) 



(4) 



and F = U,D. Here, Y^ is proportional to the quark mass matrix while Y^ gives the cou- 
plings between the heavier neutral and charged Higgs bosons and the quarks. Clearly, if Y^ 
and Y^ cannot be diagonalized simultaneously, flavor- changing neutral currents (FCNC's) 
will be induced at tree level and associated with the doublet H. If we impose some sym- 
metry to suppress the tree-level FCNC's, as in type-Il THDM, the couplings of the new 
Higgs bosons are always proportional to the quark masses. In this case, the excess Wjj 
events should be mostly b-fiavored, which is against the observation. To avoid this problem, 
instead of imposing symmetry, we find that Y^ and Y^ can be simultaneously diagonalized 
if they are related by some transformation. 

To illustrate the desired relationship between Y^ and Y^, we first introduce unitary 
matrices V^p, to diagonalize Y^ in the following bi-unitary way: 



F/a ^ vlY^VP = diag{Y^,YlY^) 



(5) 



Using 
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(6) 



where a, b and c are arbitrary complex numbers, one can easily see that 



Y, 



dia 



IfY^'^'II 



a^Y^ 
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c^Y^ 



(7) 



is still diagonal. Now if and are related by 

where /f = Vf'' IpVf and = V^^/JV^, then Y^ and Y^ can both be diagonalized by 

/F 



V£ji, as can be explicitly checked using Eq. ([5]) and the unitarity of V^m. We note that the 



matrix Ip in Eq. (E]) is not unique. More complicated examples can be found in Ref. |34|. 
Now if the quark mass hierarchy is such that Yp Yp Yp, we see in Eq. ([7]) that the 
hierarchy pattern in Y^ can be inverted with suitable choices of a, b and c. We note that since 
a, b and c are arbitrary complex numbers, all elements in Y^ are also complex in general. 
As a result, the couplings between the H doublet and light quarks are not suppressed by 
their masses. Moreover, the coupling to b quarks can be suppressed. 

To proceed the analysis, we write down the relevant interactions in terms of physical 
eigenstates: 

- 'Ch±,ho,AO = [UrT] Up + dpT] dp j -j= — 

+ (-m/?t7^Vc/l + MiVr/^rffi) + /i.e. , (9) 

where V is the Cabibbo-Kobayashi-Maskawa (CKM) matrix, and r]^ = diag(?7f , 77^^, r]^) 
contains three free parameters. For simplicity and illustration purposes, we will consider 
two schemes: 

(I): rji = rjY = vF with ^ = 1,2, 3; (10) 
(II): 7]^ = rif and r/^ = r]f for z = 1, 2, 3 and j = 1, 2 . (11) 

To suppress the coupling with the b quark, we require that 773 -C 771,2 in Scheme (I) or r/l^ <^ 1 
in Scheme (II). 

In either scheme, we search for the parameter space that can explain the excess Wjj 
events, subject to the constraint cr^jj = a{pp — )■ W H^)BR{H^ — > jj) = 4 pb, as observed 
by CDF. Moreover, we consider a 25% uncertainty in the extracted (Jwjj- Iii the scenario of 
a heavy CP-odd Higgs boson, we ignore the contribution from Fig. [11(b) and consider only 
the t-channel Feynman diagram. The contribution of Fig. [D^b) is one order less than that 
of Fig. [T^a) when ?n^o > 650 GeV. 

We first consider Scheme (I). Due to small parton distribution functions (PDF's) asso- 
ciated with charm and strange quarks in the proton (or anti-proton), we find that 772 does 



not play a significant role in determining awjj- Therefore, crwjj mainly depends on 771, the 
coupling between and quarks of the first generation, and the hadronic branching ratio of 
H^, = BR{H^ — jj). In Fig. |21 we fix 772 = 0.1. The red curves on the plane are 



contours corresponding to crwjj = (4 ± 1^ 
GeV in accordance with the CDF result 



i_pb. In this analysis, we took mass of m//± = 144 

In principle, the same t-channel diagram in Fig. [1] can contribute to Zjj events. However, 
the couplings of the Z boson to charged leptons are more suppressed than W. The blue 
curves in Fig. |2]are contours of azjj = cr{pp — )■ ZH^)BR{H^ jj) bring around 2.6 pb, 
which is the cross section of SM background process pp — t- ZZ + ZW — ?■ Zjj. We see that 
the preferred parameter region of the red curves have azjj well below the SM background. 




FIG. 2: Contours of (Jwjj = (4 ± 1) pb (thick red curves) and azjj = (2.6 ± 0.6) pb (thin blue 
curves) for Scheme (I). In this scenario, we take mf{± = 144 GeV and is sufficiently heavy. A 
K factor of 1.3 is used in computing the cross section. 

Using the extracted parameter space, we then compute the total width of using the 
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partial width formula 

r.= E ^^^H^Miivrr + ivf?] (12) 

and Sjj. Note that the 6-quark coupling has been taken to be zero in the above formula. 
When Bjj > 0.8 for 771 = 772 or 0.7 for 771 ^ 772, the total width Th± < 2 GeV, consistent 
with our narrow width approximation. This suggests that the charged Higgs boson couple 
dominantly to quarks instead of leptons. 

We now consider two cases in Scheme (II): (a) rj^ = rj ^ and (b) r]^ = O.lr]^. The 
independent parameters are then r]'^ and Bjj. Plots in Fig. |3] show that it is preferred to 
have rj^ < because it helps suppressing Zjj production. Likewise, when > 0.8, the 
total width V ii± < 2 GeV, again consistent with our narrow width approximation. 
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FIG. 3: Same as Fig. [2]but for Scheme (II). Results for rf = are shown in plot (a), and results 
for rP = 0.1?7^ are shown in plot (b). 

We note in passing that one can also consider the scenario where the roles of and 
are interchanged, with the former being heavy and the latter having a mass of 144 GeV. 
However, the parameter region for explaining the W events predicts a Zjj rate very close 
to the SM background in Scheme (I), as shown in Fig. |4](a). In Scheme (II), null deviations 
of Zjj and b-jets disfavor the small and large rf regions, respectively, as shown in Fig. IH^b). 



7 



Therefore, in comparison the previous scenario with hght charged Higgs boson and heavy 
CP-odd Higgs boson is favored. We will thus exclusively consider such a scenario in the 
following analysis of low-energy constraints. 
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FIG. 4: Contours of a{pp WA^)Br{A'^ jj) = (4 ± 1) pb (thick red curves) and a{pp — )■ 
Z^°)Br(A° ^ jj) = (2.6 ± 0.6) pb (thin blue curves) for Scheme (I) in plot (a) and Scheme (II) 
in plot (b). In this scenario, we take m^o = 144 GeV and is sufficiently heavy. A K factor of 
1.3 is used in computing the cross section. 

If the charged Higgs boson is a candidate for the new resonance, it will also induce 
interesting phenomena in low-energy systems, where the same parameters are involved. We 
find that the most interesting processes are the B — )■ Xg'j decay and the like-sign charged 
asymmetry (CA) in semileptonic Bq {q = d, s) decays. To simplify our presentation, we leave 
detailed formulas in Appendix|Xl Using the interactions in Eq. (j9]), the effective Hamiltonians 
for the 6 — )■ S7 and AB = 2 processes induced by H^, as shown in Fig. [5l are respectively 
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given by 



n 



16? (Qthiyt) + hivt)) Or, 



n{AB = 2) 



{QtJi{yt) + J2iyt))0',^ 

so 2 h{yt)qi,,PLbqYPLb 



4(47r)2m2 



(13) 



2(47r)2m^± 

where = mj/m^±, and are defined in the Appendix, = 2/3 is the top-quark 
electric charge, and 

Ji{a 
h{a 
J2{a 

z,\L — a)" \L — a)" 

1 ( \ 2 l+g , 

J^ia) = —— — — Tz —ma. 
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3 — a 


In a 


(1-a) 


2 ' (l-a)3' 


1 — 5a 


— 2a2 a^ In a 


6(1 - 


a)3 (1-a) 



1 + a ^ a In a 



2(1 -a)2 (l-a)3 

1 + a a In a 
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2 1 + a 



(l-a)2 (l-a)3 
Using the hadronic matrix elements defined by 



(14) 



{B,\qPLbqPLb\B,) 



m 



Ba 



rriBjB^Bg , 



(15) 



24 \ nib + ITT-q 

and the formulas given in the Appendix, the dispersive part of Bq-Bq mixing is found to be 

Ml^ = Mff^ + Mff = Mff^Aqe^ , (16) 



where 



1 + R% + 2Rjj cos 2^^ 
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tan0^ 
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1 + cos 2^1^ 
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FIG. 5: (a) b ^ s'j transition, where the dot indicates another possible place to attach the photon 
and (b) with q = d,s induced by H^. 

Since the charged Higgs boson is heavier than the W boson, its influence on is expected to 
be insignificant. Therefore, we set ~ Tf'.^'^^ in our analysis. Using 0^ = arg(— Mj'2/ri2)5 
the if^-mediated wrong-sign CA defined in Eq. (lASP is given by 

1 sin I ' 



a 



with 



-2/3, - 7, 



SM 
Q 



38 



and (h 



Ag sin^SM 



«!.(SM) 



(18) 



Consequently, the like-sign CA in 
Eq. dMl) is read as A^^ ^ 0.506 af^{H^) + OAMaUH^, where the SM contributions are 



38 



39| 



a^^(SM) ^ 1.9 X 10-5 and af^(SM) ^ -4.1 x 10"^ 

In the following, we numerically study the charged Higgs contributions to the B — )■ X^'j 
decay. In Scheme (I), as 773 = 77^ = r/l' <^ 1 is assumed, it is clear that their contributions to 
the B decay are small. Thus, we concentrate on the analysis of Scheme (II). Using Eqs. ( IT3ll 
and (lAip . the if ^-mediated Wilson coefficients for b s'j are given by 

M\2 



5C7 



6C' 



8V2ml±GF 



{QMvt) + Hyt)) 



mt 



{QtJi{yt) + J2{y,)) . 



(19) 



Here the enhancement mainly comes from the large r]^ coupling. Taking Eq. ( ]A3I) and 
setting 7]^ = prf^ and (pn = arg(?7^ ^7^*); one can calculate the branching ratio of -B — ?■ Xs'j 
as a function of rj^ and (pn- The 2cr range of experimental measurement B{B — )■ ^^7) = 
(3.55±0.26) X 10"^ (43| demands the two parameters to be within the shaded bands in Fig.[6l 
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where plot (a) and (b) use p = 1 and 0.5, respectively. The results show that B{B — )• Xg'j) 
is insensitive to the new phase (pn and that the allowed range of t]^ is compatible with 
the above analysis for the Wjj events. In addition, we also show in Fig. E] the constraint 
from measured ^ttlb^ (dashed blue curves). We only take into account ArriB^ here simply 
because the measurement Am^^ = 0.507±0.005 ps~^ is more precise and thus stringent than 
ArriB^ = 17.78 ± 0.12 ps~^. It is observed that the measurement of Am^^ further excludes 
some of the parameter space allowed by the B — t- Xgj decay. Finally, we superimpose 
contours of the like-sign CA (solid red curves) in Fig. |6l The like-sign CA has a strong 
dependence on the value of p. When p ~ 0(1), can be of the order of 10"'^. However, it 
drops close to the SM prediction when p ~ 0(0. 1). 

We now comment on the constraints from K-K and D-D mixings. In the usual THDM, 
contributions from box diagrams involving the charged Higgs bosons to the mass differ- 
ence are important because the charged Higgs couplings to quarks are proportional to their 



masses. Therefore, the Glashow- Iliop oulos- Maiani (GIM) mechanism (4^ is not effective to 
suppress such new physics effects 46|. In the scenarios considered in this work, the H'^qq' 
couplings are simply proportional to the CKM matrix elements. Therefore, the box diagrams 
involving the charged Higgs boson will have GIM cancellation in the approximation that the 
masses of quarks in the first two generations are negligible. Although the third generation 
fermions do not have GIM cancellation, the associated CKM matrix elements are much sup- 
pressed. In addition, the new effective operators thus induced will be further suppressed by 
powers of niw /TnH±- For example, with m/^i = 140 GeV, p = 1, rj^ = 0.4 and the dispersive 
part of K-K mixing given in Eq. (1A12I) . we obtain Amx ~ 1.58 x 10"^^ GeV, which is two or- 
ders of magnitude smaller than the current measurement, (Amji-)°^P = (3.483±0.006) x 10~^^ 



GeV 



35|. We note that unlike the conventional THDM, where the diagrams with one W- 



and one in the loop are important 45j , in Scheme (II) of our model the GIM mechanism 
is very effective in the massless limit of the first two generations of fermions. This has to 
do with the fact that the charged Higgs couplings to these quarks are independent of quark 
masses. The contributions from diagrams with the top-quark loop are also negligible due to 
the suppression of the small CKM matrix elements (Vj^V^^)^, as in the conventional THDM 
with small tan /3. The relevant formulas for K-K mixing from these contributions are given 
by Eqs. flA12p and flA13l) . The constraint from D-D mixing is even weaker in view of current 
measurements [47 1 and the fact that new physics contributions are both GIM and doubly 
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Cabibbo suppressed. 




(a) (b) 

FIG. 6: Contours of j4^^ (in units of 10~^) on the rj^-4>H plane, where the shaded band and the 
dashed curves show the constraints from the measured B{B — )■ Xg^) and Amsj^, respectively, 
within their 2a errors. 

In summary, we have studied a scenario of the two-Higgs doublet model as a possible 
explanation for the excess Wjj events observed by the CDF Collaboration. In this scenario, 
the charged Higgs boson has a mass of about 144 GeV and decays into the dijets. We find 
that both Scheme (I) and Scheme (II) considered in this work can explain the Wjj anomaly 
while not upsetting the constraints of Zjj and b-jets being consistent with standard model 
expectations. When applying the scenario to low-energy B meson phenomena, we find that 
very little constraint can be imposed on Scheme (I) as 773 couplings to the third generation 
quarks are assumed to be negligible. Scheme (II), on the other hand, has constraints from 
the B — )■ Xg'j decay and Am^^. In particular, we find that if t]'^ for the first two generations 
is of the same order of magnitude as r]^ , it is possible to obtain A^^ ~ O{10~^). Constraints 
from K-K and D-D mixings are found to be loose primarily due to the GIM cancellation. 
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Appendix A: B ^ Xsj and like-sign CA in semileptonic Bg decays 

For the B — > Xg'-f decay, the effective Hamiltonian is 

n,^s, = -^^tlVtb {cMOj, + c;(/i)o^j (ai) 

with 

87r^ 



The branching ratio is given by 40 1 



O'j, = ^.-a,.(l-75)foF^^ (A2) 



B{B^X,^)E^>i.eGeV = [aoo + a7j{\SC7\^ + \SC',\') 

+ aorReiSCj) + a'Q^Re{6Cj)] x 10"^ , (A3) 

where Ooo = 3.15 ± 0.23, 097 = —14.81, 077 = 16.68, and Oq^ = —0.23. The parameters 
6Cy = C^^ and 5Cy ~ C^^ stand for new physics contributions. 

To understand the like-sign CA, we start with a discussion of relevant phenomena. In 
the strong interaction eigenbasis, the Hamiltonian for unstable Bq and Bq states is 

H = M'i-i— , (A4) 

where (M'^) denotes the absorptive (dispersive) part of the Bq ^ Bq transition. Accord- 



ingly, the time-dependent wrong-sign CA in semileptonic Bq decays is defined and given [35 1 
by 

, _ T{Bq{t) ^ ^+x) - T{Bq{t) ^ rx) 



a; 



V{Bq{t) l+X) + V{Bq{t) i-X) 



r 



Iml^j, (A5) 
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Here, the assumption of ^ ^12 the Bg system has been used. Intriguingly, a^^ is 
actually not a time-dependent quan tity. The relation between the wrong and like-sign CAs 
is defined and expressed by [36|, [37 1 

r{bb i+i+x) - r{bb e-e-x) 



A" 



T{bb ^ i+£+X) + T{bb 
0.506(43)af^ + 0.494(43)a^, 



-X) 



(A6) 



Clearly, the like-sign CA is associated with the wrong-sign CA's of the Ba and Bg systems. 
Since the direct measurements of af^ and a^^ are still quite imprecise, either 6 — )■ or 6 — )■ s 
transition or both can be the source of the unexpectedly large A^f, observed experimentally. 
In order to explore new physics effects, we parameterize the transition matrix elements 

as 



r?2 



'-12 '^q^ 



(A7) 



for q = d,s, where 



M' 



<j,SM[NP] 



12 



A 



tan ( 



M- 



g,SM[NP] 
12 



1 + 



F 

5 12 



<j,SM 



■.g,SM[NP] 
12 



. SM[NP] 



,NP| 



12 



IM- 

12 



(j,SM| ' 



'12 

-L 1 



12 I 



rf sin2(^^^^-/3. 



1 + rf cos2(eNP_^^) 
Here, the SM contribution is 



tan 7 



r • / NP 
rg sin(7„ - 7, 



1 - cos(7NP - 7 



(A8) 



SM 



^f"''W „ _ f2 



VBmBj'BBgiV*VttySo{xt) 



(A9) 



with So{xt) = 0.784x° 



0.76 



Xt 



{mt/mwY ^'^d TjB ~ 0.55 being the QCD correction to So^Xt). 



The phases appearing in Eq. (IA7p are CP- violating phases. Note that 13d = (3d and f3s = —f3s. 



Using 

'12 



arg(— M;^2/ri2); the wrong-sign CA in Eq. ( ]A5P with new physics effects on F 



12 



and Mfn can be derived as 



\ sin </)g 



Af sin,^ 



SM 



aim) 



(AlO) 
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with = 2Pq--if^ and 



,g ^..^ = + 0^ — 7^. Furthermore, the mass and rate differences 
between the heavy and light B mesons are given by 
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AH 



2|r^2l cos I 



(All) 



As a comparison, we consider the new physics effect on K — K mixing due to the box 
diagram with the top quark and the charged Higgs boson in the loop. This is seen to be the 
major contribution as other diagrams involving lighter quarks are GIM suppressed in the 
massless limit or are smaller even when mass effects are taken into account. The result of 
Mi2 for the diagram with both the intermediate bosons being the charged Higgs boson is 



M- 



K,HH _ 



ts) 



12 



12(47r)2m2^ 



+ ^ytMvt 



1 



V 



rriK 



D|4 



t/|4 



h{yt 



5 Kms + nid 

2 

rriK 



(A12) 



4 \ms + 

The contribution from the diagram with one W boson and one charged Higgs boson in the 
loop is 



12v^7r2 

1 3 / rriK 
4 2 \ms + nid 



\r]^\ Ki{xt,XH) 



m 



w 



X2 



where Xt = mf/m^, xh = rn'^H±/mly and 

= '''"^irRft^Dx.H- (a- 6M» 

Hadronic effects have been included in Eqs. flA12|A13P already. 
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